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Most of the mass content of dark matter haloes is expected to be in the form of tidal debris. 
The density of debris is not constant, but rather can grow due to formation of caustics at the 
apocenters and pericenters of the orbit, or decay as a result of phase mixing. In the phase space, 
the debris assemble in a hierarchy which is truncated by the primordial temperature of dark matter. 
Understanding this phase structure can be of significant importance for the interpretation of many 
astrophysical observations and in particular dark matter detection experiments. With this purpose 
in mind, we develop a general theoretical framework to describe the hierarchical structure of the 
phase space of cold dark matter haloes. We do not make any assumption of spherical symmetry 
and/or smooth and continuous accretion. Instead, working with correlation functions in the action- 
angle space, we can fully account for the hierarchical structure (predicting a two-point correlation 
function cx AJ -1,6 in the action space), as well as the primordial discreteness of the phase space. As 
an application, we estimate the boost to the dark matter annihilation signal due to the structure of 
the phase space within virial radius: the boost due to the hierarchical tidal debris is of order unity, 
whereas the primordial discreteness of the phase structure can boost the total annihilation signal 
by up to an order of magnitude. The latter is dominated by the regions beyond 20% of the virial 
radius, and is largest for the recently formed haloes with the least degree of phase mixing. 



I. INTRODUCTION 



Cosmological N-body simulations show that dark mat- 
ter (DM) haloes which form in a ACDM Universe con- 
tain a large number of subhaloes of all sizes and masses. 
What remains outside the subhaloes are ungrouped indi- 
vidual particles whose masses set the resolution limit of 
the simulation. If the simulations where to have enough 
resolution to resolve every single subhalo then it is ex- 
pected that the smallest subhaloes would be the micro- 
haloes of about 1(T 6 M §,§,§. Does all the mass 
of a given halo reside inside the gravitationally bound 
subhaloes ? As a subhalo falls through the gravitational 
field of its host halo, it becomes tidally disrupted. A 
tidal stream extends along the orbit of the subhalo and 
can contain a large fraction of the satellite mass. There- 
fore, a significant fraction of a DM halo is expected to be 
in the form of streams and caustics. Depending on their 
length, the density of the streams can vary and is rela- 
tively not very large. However, as a stream folds back on 
itself, zones of higher density, i.e. caustics, form (see e.g. 
0, [| ||, 0, ||). In principle, these arc not true caustics 
but only smcared-out caustics due to finite DM velocity 
dispersion, however, it is convenient to refer to them sim- 
ply as DM caustics. Hereafter we shall refer to unbound 
streams and caustics jointly as tidal debris. 
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Dark matter tidal debris, so far mostly unresolved in 
cosmological N-body simulations, are expected to pop- 
ulate our own halo. Many stellar counterparts to such 
debris have been detected so far {e.g. || |l0|]) and many 
more are expected to be detected with future missions 
like GAIA. The hierarchical growth of the host halo from 
the disruption of satellite haloes reflects in a hierarchical 
structure of the phase space. The true lowest cutoff to 
this hierarchy is not set by the microhalocs but by pri- 
mordial dark matter velocity dispersion. The hierarchical 
phase structure indicates that after removing all bound 
subhaloes from a given DM halo, its phase space remains 
still unsmooth due to debris from disrupted subhaloes. 
The tidal debris are never smeared out because of con- 
servation of phase space density and volume, although 
they become less dense as they wrap around the halo. It 
is this phase structure which we study here. 

Secondary infall or self-similar accretion model pro- 
vides a solid theoretical base for the study of halo for- 
mation, and models the phase structure of DM haloes 
jnj 12 . However, since this model assumes continuous 
accretion, it cannot capture the hierarchical nature of 
halo formation. On the other hand, numerical simula- 
tions still lack enough resolution to resolve the hierarchi- 
cal phase structure, although progress is being made in 
this direction g, [H| fO). 

Here, we aim at capturing the hierarchical phase struc- 
ture of dark matter haloes and its intrinsically discrete 
nature, without resorting to any assumption of spheri- 
cal symmetry or smooth and self-similar accretion. We 
divide the structure of a dark matter halo into three cate- 
gories: (1) the primordial and intrinsically discrete phase 
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structure, formed prior to any merger or accretion and 
entirely due to the coldness of the initial condition; (2) 
the hierarchical phase structure of tidal debris from dis- 
rupted satellites, and (3) the hierarchical phase struc- 
ture of undisrupted subhaloes. We leave the studyof the 
undisrupted substructures to a companion paper [|l5| and 
in this work we only study cases (1) and (2). 

To study phase structure induced by debris from dis- 
rupted satellites, we assume that at a given level in the 
hierarchy, all structures added earlier and which lie at 
smaller scales are smooth. This sets the lowest level of 
the hierarchy at the scale determined by the velocity dis- 
persion of earliest dark matter haloes. However, this is 
not entirely correct since the earliest dark matter haloes 
themselves are not smooth and have a structure which is 
due to the coldness of the initial condition. Thus, there 
is a fundamental discreteness scale which is determined 
by primordial dark matter velocity dispersion (see Fig. 

This complicated process is studied here through cor- 
relation functions in the action-angle space where Hamil- 
tonian is only a function of the adiabatic invariants, i.e. 
the action variables. Their conjugate variables, the an- 
gle variables increase linearly in time. The action-angle 
variables arc extremely useful for studying tidal streams 
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is 



|19L |20f| . However, working with the action- 
angle variables, we are restricted to regions within the 
virial radius (with a quasi-static potential) , and hence the 
phase structures that might arise outside the virial radius 
(e.g. between the virial and the turnaround radii) can- 
not be studied in the present framework. For direct DM 
detection and cosmic-ray signal of DM annihilation, only 
the nearby phase structure plays a role and our method 
is valid (see e.g. pl[). However, for lensing experiments 
and 7-ray emission from DM annihilations, for example 
from other galaxies, the structures outside the virial ra- 
dius can be rather important (see e.g. |22||). 

We assume that the satellite orbits are integrable in the 
host DM potential (although, in [v] we remark on chaos 
and non-integrable systems). Therefore, the phase space 
distribution can be described in terms of the action-angle 
variables, {Ji,9i}, so that: 



(1) 

(2) 



where the Hamiltonian, 7i = H[J], is only a function of 
action variables, Ji and fi^s are the angular frequencies. 
Fig. (Q) shows a cartoon picture of phase mixing in the 
action-angle space, and its correspondence to the real 
space. 

The hierarchical phase structure and its fundamental 
discreteness set by primordial DM velocity dispersion are 
captured by the correlation function of the phase density. 
Since, after a long time, the distribution in the angle 
space is uniform, the phase density is only a function of 
the action variables. This can be easily seen by writing 




FIG. 1: The top horizontal panel shows the phase space of 
the merger of two dark matter haloes, each of which has its 
own hierarchy of phase structure. The times on the top panel 
refer to the crossing times. The zooming shows that each 
hierarchy contains a lower level and so on. The hierarchy is 
cut at the scale of the smallest dark matter halo that has 
been accreted to the final halo. However, the phase space 
is not smooth below this scale. Indeed, the phase space is 
intrinsically discrete due to the coldness of dark matter shown 
by the last zooming on the left. (Top panel: courtesy of 
Vlasov-Poisson simulation |E3[.) 
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FIG. 2: A one dimensional cartoon of the evolution of tidal 
streams in both phase and action-angle spaces. As structure 
wraps around the phase space, more streams cross the same 
angle coordinate, which leads to a discrete lattice-like struc- 
ture in the action space. 



the collisionless Boltzmann equation for the equilibrium 
distribution in the action-angle space: 



dt 







;df-df 



Ji 



0. 



(3) 



which, combining with Eq. (g), implies that the equilib- 
rium phase space density can only be a function of action 
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variables (and is known as the strong Jeans theorem p4j). 

This enables us to evaluate the density-density corre- 
lation function. Our results are only valid statistically 
for typical haloes and thus may not agree with results 
obtained for individual haloes in the simulations. 

The nature of DM remains a mystery. Supersymmetry 
and extra-dimensional extensions of the standard elec- 
trowcak model provide a natural candidate in the form 
of a weakly interacting and massive particle (hereafter 
WIMP). These species should fill up the galactic halo. If 
DM consists of WIMP's, they are expected to strongly 
annihilate in the dense regions of our halo and gener- 
ate in particular gamma-rays and charged cosmic rays. 
Hence, hierarchical structure of phase space can lead to 
the enhancement of DM annihilation signal J25| . We also 
evaluate the boost to the annihilation signal due to tidal 
debris and discreteness of the phase structure. We show 
that the boost due to tidal debris is of order one, whereas 
the boost from the discrete phase structure can be up to 
one order of magnitude higher. 

In Section [h] we review a few basic relations for action- 
angle variables, 
the correlation 
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phase structure due t o tidal debris and their discrete- 
we evaluate the boost on the an- 
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ness. In Section IV A, 

nihilation signal due to tidal debris. In Section IV B| . 
evaluate the boost of the annihilation signal from intrin- 
sic discreteness of the phase structure, and finally Section 
^ concludes the paper. 



II. STREAMS AND COHERENCE VOLUME OF 
THE PHASE SPACE 

We use the definition of action-angle variables (^J) and 
assume that the frequencies arc not degenerate, i.e. the 
Hessian matrix 



Hi 



d 2 n 

dJjdJi 



(4) 



has non-zero eigenvalues, or cquivalently, a non- vanishing 
determinant: 



(5) 



with the possible exception of a zero measure region of 
the phase space. Note that, this implies that the halo 
potential cannot be assumed to be exactly spherically 
symmetric, as two of the frequencies would be equal. 

If a satellite galaxy has originally a small spread in the 
action variables, AJ^, its spread in the angle variables 
increases as: 



A9i = (HijAJj) t acc , p + A9q, 



(6) 



where t acCiP is the time since accretion of the progenitor 
of the debris into the halo. The last term, the initial 
extent of the debris, is subdominant at large times. We 
set this term to zero for now, but at the end of Section 



IV B, wc discuss when it can become important and how 



it could affect our results. Therefore, the total volume 
swept in the angle space grows as 



A 3 9 = (|^,|A 3 J) tl cc , p = (A 3 fi) t 3 QCC , p 



(7) 



where we used the definition of TLij in equation (Q), and 
A 3 f2 is the volume occupied by the debris of satellite 
particles in the frequency space. 

As the total volume of the angle space is (27r) 3 , the 
number of streams passing through each angular coordi- 
nate is 



iVstre 



(A 3 ft) 



2- 



(8) 



Thus, the total mass of each stream, m s t r eam, is the mass 
of the debris, to, divided by A^ stream 



\A 3 nJ 



t a 



cc,p 



2 71 



(9) 



Put another way, the action space is divided into cells of 
volume: 



A 3 J st 



2tt 



(10) 



as a result of phase space mixing {e.g. Q) . 

With this picture in mind, we can write the distribu- 
tion in the action space as the sum of the contributions 
from individual progenitors: 



f(J,0) — y^/p(J,fl;£q C c,p), 



(11) 



where each f p has a cellular structure characterized by 
Eq. (|l0|), as shown in Fig. (^), which gets finer and finer 
with time. Eq. ( JTT| ) is the phase space analog of the 
widely used halo model in cosmology [^6|, [27j , where the 
density is assumed to be the sum of contributions from 
individual haloes with given profiles. Correspondingly, 
f p characterizes the profile of individual progenitors in 
our picture. 

We can now write the real space density as: 

p(x,t)=Y / f d 3 Jd 3 8f p (J,9;t acc , p )6 3 D [x-x(J,9)} 
v J 

= Y1 [ d 3 J/ p (J,x,i acc , p )p(x;J),(12) 
p J 

where /5(x; J) is the density of a distribution of unit mass, 
with a fixed action variable J, and uniform angle distri- 
bution: 



p(x;J) 



d 3 e 

(27T) 



3 ^[X-X(M)], 



(13) 



while 



/p(J,x,t acCiP ) = j5(x; J)" 



d A 6 



/p(J,0,W)<% [x-x(M)] 
(14) 
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FIG. 3: The action-space distribution of debris in a unit 2d torus with unit particle mass and no potential. The debris is 
originally within < x, y < 0.1, and —10 < v x , v y < 10. The figures show a cut through the action space with 0.09 < x,y < 0.1, 
which is characterized by / P (J, x, t aC c, P ) (Eq. 0) in our formalism. 



An example of f p is shown in Fig. (J3|) for debris in a toy 
model of a unit torus. As we will explicitly show in IV A , 
projecting this discrete structure in the action space of 
the debris into the real space leads to discrete, (nearly) 
singular, caustics that are only smoothed by the original 
velocity dispersion of the progenitor. 



III. CLUSTERING IN THE PHASE SPACE 

Averaging over different possible realizations of the de- 
bris within a halo, the mean phase space density can be 
written as an integral: 

(^2fp(J, x ,tacc,p)) = J dAp.g (1) (J-J p ,x,t acC:P ), (15) 



where and 3 P are the profile and mean action of 
individual progenitors, while 



dn 



dN p _ dm p d 3 J p dt acc , p d , 3j dt 



(16) 



pu^acc^p 



is the differential progenitor number density per units of 
progenitor mass, m p , it action space volume d 3 J p , and its 
accretion time t acCtP . We now follow an analogy with the 
cosmological halo model pi, p7| to write the clustering 



in the action space as a superposition of one and two- 
progenitor terms: 

( ^ t fpi (Jl i x j tacc,p 1 )fp 2 (J2, x, t acCtP2 ) \ 

\Pl,P2 / 



dN p {l - prog.) + / dN pl / dN P2 {2 - prog.). (17) 



The one-progenitor term characterizes the self-clustering 
of individual progenitor action-space profiles. 



(1 - prog.) = g ( >(J 1 -J p ,TC,t acc . p )g ( > (J 2 - J p , x, t acc . p ), 

(18) 

while the two-progenitor terms characterize the correla- 
tion between phase space density at different action vari- 
ables, within different progenitors: 



(2- prog.) = g {2) (3i - J Pl , J 2 - J P2 ,x,i 

-,(2) 



! tacc.pi ) 



— Scon. (J 1 Jpi i J2 Jp 2 1 x i ^acc,px , t a cc,p 2 ) 

(19) 

In the limit that the mean actions of different progenitors 
are not correlated, the connected part of the (2-prog.) 

term goes to zero: g C on. 0, and thus the two-progenitor 
term reduces to the correlation within the smooth halo. 
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Note that this limit cannot be strictly realized, as due to 
phase space conservation, phase streams tend to avoid 
each other, leading to < at small separations 

J Pl — J P2 . However, Liouville's theorem is not valid 
for coarse-grained phase space density, and thus coarse- 
grained progenitors can overlap in the action space. 

The connected part of the two-progenitor term origi- 
nates from the clustering of the initial conditions of the 
progenitors of the host halo, which is generally expected 
from the clustering of cosmological haloes. However, the 
structure of the one-progenitor term is more subtle: In 
addition to the cellular structure described in the previ- 
ous section (Fig. ||), the internal structure of each pro- 
genitor prior to its accretion onto the host halo would 
introduce a hierarchy within each cell. In fact, in a hier- 
archical picture of structure formation, one expects the 
sub-cellular structure of the one-progenitor term to be in- 
herited from the two-progenitor terms within progenitors 
prior to their accretion onto the main halo (see Fig. [l]). 
The key difference between the two hierarchies, however, 
is that phase mixing only continues in the action space 
of the main halo, and (following the tidal disruption) has 
stopped in the action spaces of the progenitors. 

For statistically self-similar initial conditions, we ex- 
pect the sub-cellular and two-progenitor terms to blend 
into one roughly self-similar structure in the action space, 
although individual realizations have periodic structures 
with the characteristic volume given in E n. (|lC ) . We pro- 
vide a scaling ansatz for this structure in III A. However, 
the self-similarity is cut-off by the free streaming of dark 
matter particles on small separations, due to their finite 
intrinsic velocity dispersion. This is responsible for the 
fundamental discreteness of the phase space distribution 
(see Fig. [l]), which we model in [II B . 



A. hierarchical phase structure from tidal debris 

We first consider the phase structures due to tidal de- 
bris. Once again, we emphasize that these are the the 
tidal streams that have fallen into the gravitational field 
of the host halo and are no longer bound to the original 
satellite. 

The first level of approximation that we will use to 
study phase space clustering of cold dark matter (CDM) 
is to assume a (statistically) hierarchical formation his- 
tory, where any trace of the cold initial conditions has 
been wiped out through phase mixing. Furthermore, we 
ignore the possibility of gravitationally bound structures 
in this paper (see the companion paper (lq| on this sub- 
ject). The impact of cold initial conditions will be ad- 
dressed in subsequent sections. 

Assuming uniform distribution in angles (or complete 
phase mixing), the density at each point in the halo is 
given by: 



(20) 



where /(J) is the phase space density, while p(x; J) was 
defined in Eq. (|l3|) . Note that in ( |20| ) , the function p has 
the dimension of inverse volume, 1/V(J). 

We shall assume adiabatic invariance; the action re- 
mains constant as new structures are added on larger 
scales. Hence, the distribution function in the action 
space, /(J), does not change with time, except when new 
structures are added due to satellites that are newly ac- 
creted inside the virial radius. 

The assumption of uniformity in the angle space al- 
lows us to separate the effect of phase mixing from that 
of hierarchical structure formation. While the former is 
the cause of original caustic formation, too much phase 
mixing (within a fixed potential) will eventually smooth 
out the real space density distribution. The assumption 
of a smooth /(J) distribution, implies that phase mixing 
is complete. 

On the other hand, the effect of hierarchical struc- 
ture formation is captured in /(J), through the fact that 
structure in /(J) is added on different scales, at different 
times in the history of the halo. A statistical measure 
of this history is the two point correlation function of 
the action space density. We thus hypothesize that the 
correlation function: 



0(Ji,J 2 )^(/(Ji)/(J 2 )) 



(21) 



should be a power law for statistically self-similar initial 
conditions: 



(/(Jl)/(J 2 ))dcbris^A|J|- Q |Ji- J 2 



(22) 



as long as |Ji - J 2 | < |J| with J = (Ji + J 2 )/2. The 
form of the correlation function uses the symmetry be- 
tween Ji and J 2 . It also guarantees that small scale 
structures are captured, as structures are added on dif- 
ferent scales at different times (see Figs. |l| for demon- 
stration). Moreover, since actions remain constant in the 
adiabatic invariance approximation dA/dt ~ on small 
scales. In other words, the correlation function £/(Ji, J 2 ) 
grows inside-out in the Ji, J 2 space. 

The mean phase space density of a virializcd halo, as- 
suming a virial overdensity of ~ 200, is given by: 

(2^ 3 J~i4~5 (23) 



(/(J)) 



10H 
10H 

G 2 M' 



(24) 



where M is the halo virial mass. The virial action vari- 
able is also roughly: 

Given that M cx a 6 /("rff+ 3 ) an d H cx a" 3 / 2 when per- 
turbations grow during the matter-dominated era, with 
a being the cosmological scale factor, and n e g the slope 
of the linear power spectrum, we conclude: 



3("cff+ 7 > 

(/(J)) OC J ^TTT 



P (x) = (2^ / <fJ/(J)p(x;J), 



3Kff + 7) 
n e e + 11 



1.6 ±0.1, 

(26) 
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where we have assumed n e ff ~ —2. 5 ±0.5 for cosmological 
haloes. 

An alternative way to derive ( p6| ) is to consider the 
self-similar collapse models of Fillmore & Goldreich , 
where they calculate the actions and use adiabatic in- 
variance to find the outcome of spherical cold secondary 
infall in an Einstein-de Sitter universe. For the spherical 
self-similar linear initial condition of 



(27) 



they find the action at the turn-around radius scales as 



J ta ocM t £ + 1/3 **-i,M ta oc** 



(28) 



for e < 2/3, where M ta is the mass within the turn- 
around radius. Although ( p?8| ) is only for the radial ac- 
tion, and the two other action variables vanish due to 
spherical symmetry, one may imagine that for triaxial 
CDM haloes, the three actions would become compara- 
ble: Jtf, ~ Jg r~j eJ r r~j eJ ta , where e characterizes the 
triaxiality of the halo (not to be confused with the self- 
similar profile index e in (p7j)). Therefore, eliminating 
time from the two equations in (pq), the phase space 
density, /(J) scales as: 



/(J) 



Mtg 
J 3 



(29) 



Assuming that the self-similar linear density profile has 
the same radial/mass scaling as the variance of the cos- 
mological density fluctuations, a(M) cx M _ ^ eff+3 ^ 6 
yields e ~ (n cff + 3)/6 (< 2/3 for CDM Harrison- 
Zcl'dovich primordial power spectrum), which, plugging 
into (|29|), reproduces (26|). 



B. fundamental discreteness of the phase space 
structure 



In the previous section, we considered the hierarchical 
addition of tidal debris to a DM halo. However, the hi- 
erarchy has a lower cut-off set by the velocity dispersion 
of the smallest accreted satellite. Micro haloes of ~ 10 -6 
solar mass could indeed determine such a cut-off ||, [| . 
However, this cutoff is far above the primordial velocity 
dispersion of DM itself. Therefore, the primordial veloc- 
ity dispersion introduces a fundamental discreteness in 
the phase structure. In other words, the smooth phase 
space distribution of the last section ignores the discrete 
nature of multiple streams in the phase space due to the 
presence of a cut-off in the CDM hierarchy. This dis- 
creteness shows up as a cellular or lattice structure in 
the action space, with a characteristic cell volume given 
in ( [To|) (see Figs. || and || for Id and 2d cartoons; More 
realistic simulated examples are discussed in |[9|). After 
averaging over different spacings, expected for different 
accretion times of different debris, the discreteness would 



only show up as the zero-lag of the action space correla- 
tion function: 



(/(Ji)/(Ja))dta 



Si 



Ji - J 2 ) (30) 



where TO s t roa m and A 3 J s t r0 am were defined in equations 
(gjuj), and here we have assumed a zero initial tempera- 
ture for CDM particles. A finite CDM temperature will 
smoothen the delta function, as the phase space density 
cannot exceed its primordial value (see Fig. [I] for a car- 
toon). 



IV. EXAMPLE: DARK MATTER 
ANNIHILATION MEASURE 

A. DM annihilation in tidal debris 

In this subsection, we consider the enhancement in the 
expectation value of the annihilation measure due to hi- 
erarchical structures built in the phase space from tidal 
debris. 

For a uniform distribution in the angle space, the ex- 
pectation value of the annihilation measure is given by: 



(P(x) 



= (2tt) 6 J d 3 J 1 d 3 M f (3 1 ,3 2 ) J d 3 xp(x;J 1 )p(x;J 2 ). 

(31) 

We remark that the above integral can also be relevant for 
the direct detection of DM, as it quantifies the variance 
of the density field. 

In order to investigate the impact of caustics, near the 
apocenters and pericenters of the orbits, we make a sim- 
ple analogy with a one-dimensional harmonic oscillator: 



H 



1 



(pi + p 2 v +pI) + 



2 2 
U) X , 



(32) 



For concreteness, we also assume the other two dimen- 
sions are compact with the length L y and £ z , although 
the Hamiltonian has no explicit dependence on y and 
z coordinates. As the evolution in the three spatial di- 
rections decouple, we can simply read off three action 
variables from the areas of phase diagrams for each di- 
rection: 



Px 1 2 



Jy — 



LyPy 

2tt 



Jz = 



LzPz 

2tt 



(33) 
(34) 



From these relations, we can find p(x; J) using its defini- 
tion in Eq. (||): 



p(x; J) 



X)' 



(35) 



7 



where 



V LyLz (x max #min)- 



(36) 
(37) 



Notice that the square root singularity in the projection 
kernel p(x; J) is very similar to the singularity expected 
near CDM caustics. However, for a smooth distribution 
in the action-space /(J), the real space density p(x) is 
an integral over the kernel (Eq. |2(J), which would lead to 
a smooth p(~x). Therefore, a discrete distribution in the 
action space is necessary to produce caustic singularities 
in the real space (otherwise known as fold catastrophes 
or Zel'dovich pancakes). 

We then notice that the toy model of Eq. ( (32| ) is similar 
to the motion in a nearly spherical potential, in the sense 
that the motion in one direction (x or radial) is limited by 
requiring constant action variables, while the two other 
directions (y and z, or angular directions) are compact. 
Based on this analogy, we will use: 



p(x; J) 



t (J)-r min (J)] 



V(3)y/[r - r min (J)][r max (J) - r] 



(38) 



where we have assumed an integrable nearly-spherical po- 
tential, with small angular momentum (and third inte- 
gral), while V(3) is the spatial volume occupied by the 
stream of action J. The radii, r m j n and r max are the 
minimum and maximum radii of all orbits with the same 
action variable. However, the general structure of the 
singularity close to boundaries does not change in other 
geometries. 

Using Eq. (|3|), we can evaluate the x integral in the 
emission measure (pi]). For Ji ~ J 2 , the integral is loga- 
rithmically divergent around r ~ fmax(Ji) — ?"max(J2) a s 
well as r ~ T'min(Ji) — ^min(J2)- Focussing on the outer 



caustic r„ 



we find: 



d xp(x; Ji)/5(x; J 2 ) 



xcosh 



c(Ji) + r m 



c(J 2 ) 



This yields: 



d d a:p(x;J 1 )p(x;J a )~ V(J)" 



c(Jl) 



hi 



c(J 2 ) 



F'AJ,; 



(39) 



, (40) 



where F 1 ~ i|J|<91nr max /9Ji and can be calculated, given 
the gravitational potential of the host halo. Therefore, 
the annihilation measure takes the form: 



$ ~ A / d 3 J|J| 



d 3 A3\A3\ 



In 



FIAJI 



(41) 



We see that since 3 — a > 0, the integral is finite and 
dominated by large A J's. The boost to the annihilation 
signal, which is introduced by small scale clustering in 



the action space (i.e. that a > 0) is thus given by: 



B, 



debris 



^smooth 



/ d 3 3 V{3)" 1 3 a J d 3 A3 AJ 


<:v 


In 




AJ\\ 


Jd 3 3V{3)-' L 3- 2a fd 3 AJ 


In 




AJ\\ 
J l J 
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(a -3) 2 



F a = 0(1), 



(42) 



given that F ~ d In r max /<9 In J is a dimensionless number 
of order unity. 

Therefore, we see that the boost factor obtained here 
for the tidal debris is dominated by larger separations, 
as seen from expression (|42] ) whereas our approximation 
( p8| ) is valid for small separations. To emphasize, (|3£) is 
valid in the vicinity of caustics whereas the integral (42j) 
demonstrates that most contributions come from large 
separations in the action space. 

To summarize, while we predict an 0(1) boost in an- 
nihilation signal due to (finite separation) clustering in 
action space, the main effect comes from large structures 
(and not caustics) which are not accurately captured in 
our framework. In the next section, we will address the 
impact of the discreteness of the phase space of CDM 
haloes. 



B. DM annihilation boost due to discreteness of 
phase structures and catastrophes 

The primordial velocity dispersion of DM induces a 
fundamental discreteness in the hierarchical phase struc- 
ture and can enhance the annihilation signal. The emis- 
sion measure due to this discreteness is calculated by in- 
serting ( |30|) in expression (|3l|). Summing over all streams 
and subtracting the smooth part which is obtained by 
smoothing over the fundamental streams, we obtain 



E 



J d X [p(x; Jstrcam) P( x j ^stream ) S mooth] 



K 



d d n 



dM, 



bale 



d 3 n 



J d 3 x [p(x; J) 2 - p(x; J)Looth] - 



(43) 



for the enhancement of the emission measure due to dis- 
creteness of the phase structure, where /5(x; J) S mooth is 
the stream density, smoothed to the level that different 
streams overlap, and for the stream mass we assume: 



^stream ^ ^ 



dAh 



halo 



d 3 n 



(44) 
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where d 3 Q is the volume element in the frequency space 
and factor K is the ratio of the density of debris to that 
of the host halo: p d cbris = m/v — K p ha io- 

The above expression is to be compared to the annihi- 
lation measure for the smooth halo: 



smooth 



,3 / rfMhalo 



3o' ( <^halo \ 



d 3 n 



\ d 3 n> J 



J d 3 xp(x;J)p(x;J'), 



(45) 



where l~i and fl' are functions of J and J' respectively. 
Therefore, the boost associated with a given point in the 
action space is given by: 



^smooth 
-3 I AMhak, 

K 2 { — \ ^ d3n 
2ir 



Jd 3 x [p(x;J) 2 -p(x;J)L ooth ] 
Jd^'(^)JcPxp(x;3)p(x;J') 



(46) 

In order to estimate the boost factor, we should first 
approximate the density integral J d 3 xp(x;J)p(x;J'). As 
we discussed in the previous section for |J — J'| <C |J| the 
integral is dominated by the regions around the turn- 
around radii (or caustics) and thus grows as 



d 3 xp(x; J)/5(x; J') 



In (I J - J' 



V(J) 



(47) 



as seen in Eq. (p0|). In the opposite limit, |J — J'| ~ |J|, 
assuming that the two density kernels overlap, p(x, J) 
can be approximated as roughly constant, which yields: 



J d 3 xp(x; J)p(x; J') ~ V _1 (J). 



(48) 



Now we note that in (^6|) , the integral over fl' (or equiv- 
alently J') in the denominator is dominated by the large 
values of |J — J'|, while the numerator is in the small 
jJ — J'| regime. Therefore, substituting from ( (47j - f48| ) , we 
find: 



B dis [J]~K 2 



-K 2 



t 

1 ( t 



d In Mi ia io 



3 V27T 



d 3 n 

dlnM, 



halo 



d 3 n 



in 



in 



IAJ 



CDM| 



A J int— stream | 
/CDM 



halo/ 



(49) 



where |AJcdm| and | AJi„t_ s troam| characterise the fun- 
damental CDM stream width and the inter-stream spac- 
ing, respectively. The stream thickness and spacing in 
real space have been calculated within the framework of 
selfsimilar model |2£|. However to obtain the ratio we 
use a far simpler approximation. To get the last line of 
expression ([|9|), we have used the fact that the volume 
in phase space occupied by the fundamental streams is 



(27t) 3 A 3 Jcdm = Mhaio//cDM, while the total volume in 
phase space of the entire halo is: (27r) 3 A 3 Ji nt _ s troam = 

Mhalo/ (/halo)- 

The logarithmic enhancement factor is roughly: 



In 



/cDM 



hale 



20 + In 



100 GeV 



(50) 



for CDM particles of mass m x . In arriving at ( |50| ) 
we have used that / C dm ~ ^Wcrit (T x /T C mb)V°cdm 
and (/halo) ~ 10 3 fi m/ 9 cr it/cr 3 ir , with T CM b ~ 10~ 4 eV, 
T x = ■|"i x (Tq DM ~ m x /40 is the CDM kinetic decou- 
pling temperature. 

In order to estimate the boost factor in (fi"9"|), we also 
need to find the volume element occupied in the fre- 
quency space, d 3 f2, by a given mass element. First we 
should notice that for a Keplerian potential ip(r) ex —r , 
this volume element vanishes as three frequencies are 
equal, i.e. £l r = fi^ = f^, where SI3 is the frequency 
for the third integral. 

For a general potential, ip(r), we have 



— - 1 = 3 + dln — - 1 
{Id, V d\nr 



(51) 



for a nearly circular orbit and where ' indicates first 
derivative w.r.t. r, i.e. ip' = dip/ dr. Moreover, the 
extent in f2 3 depends on the triaxiality of the halo: 



dQ 3 ~ ef^, 



(52) 



where e ~ 10% characterizes the typical triaxiality of 
CDM haloes. Therefore, we will approximate the volume 
element as 



where 



d 3 n ~ 4tt en 




GM h 



alo 



r 



1 dfy, (53) 



(54) 



Hence, for a general potential we obtain 
— cPti — =-^( 2 + dlnip /dlnr)(ip /r) 



y/3 + d\ntp'/d\-a.r - 1 ) {ip' /r - ip") 



(55) 

To evaluate ( (49| ) we also need to know the multiplication 
factor K, which is the ratio of the density of disrupted 
satellite to that of the host halo. This factor can be 
evaluated for a general spherical potential, assuming that 
the debris start with the maximum density that allows 
them to be tidally unbound: 



K ~ 3 {ip 1 - rip") (2(p' + rip") 1 



(56) 
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Putting (B5h and (B6|) in (|9J) wc obtain 



B, 



(lis 



36tt 2 (v9' - rp")(V + VT 1 (^3 + <^ - l) 



x (^Ar 3/2 



87rGp c 



-3/2 



— In ^ CD: 



3e V (/halo) 



(57) 



where </3 = e?(ln<//)/<i(lnr) and p cv n = |^ is the critical 
density of the Universe. Moreover, we have used the fact 
that the product of the present-day Hubble constant and 
the age of the Universe is unity (Hot = 1.03 ± 0.04), for 
the current concordance cosmology. 
For a power-law potential: 



(p(r) 



-<p r 



0*1, 



(58) 



one can see that (p7t) yields 



72tt 2 Q3 + 2)^2(1-/3) / / CDM V 3 Phato \ ~ 3/2 
Sdis = — ; — - — m( 



3c 



(/hal o ) / V Pcrit 



(59) 

where phaio, the local halo density for a general potential 
is 



P(r) = ^ + 2<J/r) 



(60) 



from Poisson equation. 

For the power-law potential ( pq ) , we can find the local 
density in terms of the critical density and the radius and 
consequently obtain the boost factor ( |57j ) as a function 
of r/r v j r where r V i r is the virial radius of the halo, defined 
as the radius within which the mean density is 200 / o cr it. 
The local boost factor for a power-law potential, and our 
nominal values of e and /cDM/(/haio), is 



36 2 

Brik = — tH 



(3 + 2 



5 (1-^)^2^5-1 Vr 



3(/3+2)/2 



(61) 



The above expression shows that the boost is most sig- 
nificant in the outskirts of the halo, and for large values 
of 0. Indeed, the local boost diverges as /3 — > 1. We note 
that, in the context of an NFW potential p9[: 



(^NFW 



-tpo- 



ln(l + x) 



(62) 



i.e. (3 ~ 1 corresponds to the outskirts of the halo and 
/3 ~ — 1 to the central part. 

Similarly to the power-law potential, wc can obtain the 
local boost factor @ for an NFW potential (|§). The 
boost is shown in Fig. [I] as a function of r/r s and for dif- 
ferent values of the concentration parameter c = r v i T /r s . 
The boost increases as we go towards the outskirts of the 
halo as the number of streams decreases, hence increas- 
ing the density of individual caustics. As we increase the 
concentration, the central density of the halo becomes 
large, which in turn decreases the local boost factor. 



Having evaluated the local boost, we can evaluate the 
total boost from the halo which wc defined as 



where 



and 



D _ 'J'total -. 

-Dtotal = ^ J- j 

^smooth 



^smooth = Tg p 2 s / 4:TTX 2 p(x) 2 d,X, 
JO 



(63) 



(64) 



*totai = ript [1 + B dis (x)] p(xY4TTx'dx, (65) 
Jo 

where B^is is given by (^) and the density profile of the 
smooth halo is given by p(x). For NFW density pro- 
file p = p s /[x(l + x) 2 ] where p s is the scale density, we 
have plotted the variation of the total boost with the 
concentration parameter, c, in Fig. (^). Again, the boost 
decreases as we increase the concentration, since the flux 
becomes dominated by the central part of the halo. 

We saw that the boost is mostly in the outskirts of 
the haloes, namely beyond the 20% of the virial radius. 
However, one has to be cautious, since in the outskirts 
of the halo the gravitational field of the halo approaches 
a Kcplcrian potential, causing the frequencies to become 
degenerate and the term A#o in which wc have so 
far ignored, can become important. Thus, we need to 
study the importance of this term for our analysis and 
the boost. Expression (0) now becomes 



1 + -^ 

t An 



(66) 



Hence the volume element in the angle space (0) should 
be replaced by 



A 3 (9 = (A 3 n) t 3 det 6.. 



ld9 



t dVL 3 



(67) 



where det stands for the determinant. For circular orbit 
approximation we have tt^ = y ip' /r and we use expres- 
sions (^lj) and (|5^) for volume element in the frequency 
space. We thus find that the boost ( |57|) has to be multi- 
plied by the inverse of the following determinant: 

det(*« + H?£)-ll 
x [ 1 + 




if the frequencies were to become near degenerate. In 
obtaining (^) we have used A#o / AO ~ 1/Q . The effect 
of this factor in reducing the boost is shown for an NFW 
potential (^) by the dashed lines in Figs. ^ and |j| The 
local boost is reduced slightly in the outskirts as expected 
and the total boost is reduced by a factor of about 2. 
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FIG. 4: The local boost factor due to primordial discreteness of the phase structure, for an NFW potential: The plots show how 
the boost in the annihilation measure of a DM halo changes as we go from the inner part of halo to outer parts, due to the 
discrete phase space structure of CDM. The local boost increases as we go towards the outskirts of the halo and also as we 
decrease the concentration. The dashed curves on the left panel show the boost if we include corrections due to a finite initial 
phase for the debris (pq) , which become important for nearly degenerate frequencies. The right panel shows that most of the 
boost comes from regions beyond 20% of the virial radius. 




concentration 

FIG. 5: The estimated total boost in the annihilation measure 
of a DM halo, due to the discrete distribution in the CDM 
phase space is shown for an NFW halo. The lower dashed 
curve shows the total boost when we include corrections due 
to a finite initial phase for the debris (^), which become 
important for nearly degenerate frequencies. 



V. CONCLUSIONS 

Working in phase space and with action-angle vari- 
ables, we have shown that the density-density correlation 
function can capture the hierarchical phase structure of 
tidal debris and also the fundamental discreteness of the 
phase structure due to the coldness of the CDM initial 
conditions. The study presented here assumes no spheri- 
cal symmetry, no continuous or smooth accretion, and no 
self-similar infall for the formation of dark matter haloes. 
It is thus a general scheme for quantifying the statistical 
properties of the phase structure of the virialized region 
of cosmological haloes. 

As an application, we have obtained the significance 
for dark matter annihilation signal due to the hierarchi- 
cal phase structure of tidal debris and have shown that 
this structure boosts the annihilation flux by order unity. 
On the other hand, the total boost due to the primordial 
discreteness of the phase space can be one order of mag- 
nitude higher for low-concentration (or recently formed) 
dark matter haloes. 

While this paper dealt with unbound debris and caus- 
tics in dark matter haloes, in a companion paper JlS} ], we 
calculate the boost to the annihilation signal due to the 
gravitationally bound substructure or sub-haloes. Com- 
bining the results of both papers, we can write down a 
concise and approximate formula for the local boost due 
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halo 



FIG. 6: The estimated local total boost including contribu- 
tions from the debris, discreteness and the subhaloes, given 
by ([ji]) with the first term set at its lowest value of unity, is 
shown for different redshifts. At high redshifts, the primordial 
caustics dominate over all other effects. However, at low red- 
shifts the discreteness effect due to caustics is only important 
in the outskirts of the haloes. 



to all substructures: 
(p 2 (x)) 

Boost = , , „ — 1 = -Bdobris + £>dis + -B su b 



(P(x)) 2 
0(1) + 3x10 



5 / Pcrit 



3/ + 10 6 (_***)( * 



(69) 



which should be valid within a factor of 3 in the virialized 
region of the haloes, phaio is the local coarse-grained 
density of the halo at redshift z, while p CI - lt is the critical 
density of the Universe at redshift z. The first term in 
( |69| ) is due to the hierarc hical s tructure of CDM debris, 
that we estimated in Sec. IV A. The second term is due 



to the discrete nature of the CDM phase space, where we 
used ([59]) with (3 ~ 0.5 and our nominal values for other 
parameters. Lastly, the third term is the contribution 
due to gravitationally bound sub-haloes pa]. We take 



H 2 /Hg = VL m {l + zf + fi A and plot (f|) in Fig. §. 

One may wonder whether the discreteness of the phase 
space of sub-halos could lead to an additional boost in the 
annihilation signal. In other words, should we add -B su b 
and £?dis to get the total boost, or rather should they be 
multiplied? To answer this question, we notice that the 
main contribution to -B S ub is due to the smallest sub n - 
halocs (or micro-haloes) which have the highest densities 
|fl5f , while the B^is is mainly due to the lowest density re- 
gions of the haloes, which have the lowest degree of phase 
mixing. Therefore, we expect the two terms -B S ub and 
£?dis to simply add incoherently, as the cross-correlation 
between the two sources of sub-structure should be small. 

Finally, we should point out that the results here apply 
to phase structure within virial radius, and those outside 
the virial radius which we have not studied here, might 
yield a bigger boost factor. It is reasonable to study 
the streams and caustics that lie between the virial and 
the turnaround radii by using the secondary infall model 
p| , as radial approximation can be reasonably applied 
to regions beyond the virial radius (see Fig.[|). 

Finally, we remark on the most instrumental assump- 
tion in our framework, which was the integrability of or- 
bits in the CDM potential, since one cannot define action- 
angle variables in a non-integrable system. This is char- 
acterized by the appearance of chaotic orbits in parts of 
the phase space. First, we should point out that as CDM 
haloes have a triaxial structure, a significant fraction of 
halo particles cannot be on chaotic orbits. Moreover, 
the difference between chaotic and integrablc orbits only 
becomes important after many orbital times, which are 
only possible in the inner parts of the halo. As most of 
the boost to the annihilation caused by the discreteness 
in the action space comes from the outskirts of the halo 
(Fig.^), we do not expect a significant difference due to 
chaotic orbits. Nevertheless, the implications of chaos for 
the structure of the CDM phase space correlation func- 
tion remains an intriguing question. 
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